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I— H ■ Abstract: General high-order rogue waves in the nonlinear Schrodinger equation are derived by the bilinear 

' method. These rogue waves are given in terms of determinants whose matrix elements have simple algebraic expres- 

c*j . sions. It is shown that the general iV-th order rogue waves contain N — 1 free irreducible complex parameters. In 

* *~| ' addition, the specific rogue waves obtained by Akhmediev et al. (Phys. Rev. E 80, 026601 (2009)) correspond to 

Ch . special choices of these free parameters, and they have the highest peak amplitudes among all rogue waves of the same 
order. If other values of these free parameters are taken, however, these general rogue waves can exhibit other solu- 

y—i . tion dynamics such as arrays of fundamental rogue waves arising at different times and spatial positions and forming 

' interesting patterns. 
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^ • 1 Introduction 
O 



Rogue waves, also known as freak waves, monster waves, killer waves, extreme waves, and abnormal waves, is a hot 
topic in physics these days. This name comes originally from oceanography, and it refers to large and spontaneous 
ocean surface waves that occur in the sea and are a threat even to large ships and ocean liners. Recently, an optical 
analogue of rogue waves — optical rogue waves, was observed in optical fibres [1] [2]. These optical rogue waves are 
narrow pulses which emerge from initially weakly modulated continuous-wave signals. A growing consensus is that 
both oceanic and optical rogue waves appear as a result of modulation instability of monochromatic nonlinear waves. 
Mathematically, the simplest and most universal model for the description of modulation instability and subsequent 
nonlinear evolution of quasi- monochromatic waves is the focusing nonlinear Schrodinger (NLS) equation [3J El E] ■ This 
equation is integrable [6], thus its solutions often admit analytical expressions. For rogue waves, the simplest (lowest- 
order) analytical solution was obtained by Peregrine [7j. This solution approaches a non-zero constant background 
as time goes to ±oo but rises to a peak amplitude of three times the background in the intermediate time. Special 
higher-order rogue waves were obtained by Akhmediev, ct al. using Darboux transformation [8]. These rogue waves 
could reach higher peak amplitude from a constant background. Recently, more general higher-order (multi-Peregrine) 
rogue waves were obtained by Dubard, et al. [9j [10] and Gaillard jTT] . It was shown that these higher-order waves 
could possess multiple intensity peaks at different points of the space-time plane. These exact rogue-wave solutions, 
which sit on non-zero constant background, are very different from the familiar soliton and multi-soliton solutions 
which sit on the zero background. These rogue waves are intimately related to homoclinic solutions [T2]. Indeed, 
rogue waves can be obtained from homoclinic solutions when the spatial period of homoclinic solutions goes to infinity 



*Email: ohta@math.kobe-u.ac.jp 
t Email: jyang@math.uvm.edu 



1 



[T3"] . These rogue waves are also related to breather solutions which move on a non-zero constant background with 
profiles changing with time [14]. 

In this article, we derive general high-order rogue waves in the nonlinear Schrodinger equation and explore their 
new solution dynamics. Our derivation is based on the bilinear method in the soliton theory |15j . Our solution is 
given in terms of Gram determinants and then further simplified so that the elements in the determinant matrices 
have simple algebraic expressions. Compared to the high-order rogue waves presented in [HUH], our solution appears 
to be more explicit and more easily yielding specific expressions for rogue waves of any given order. We also show 
that these general rogue waves of 7V-th order contain N — 1 free irreducible complex parameters. In addition, the 
specific rogue waves obtained in [5] correspond to special choices of these free parameters, and they have the highest 
peak amplitudes among all rogue waves of the same order. If other values of these free parameters are taken, however, 
these general rogue waves can exhibit other solution dynamics such as arrays of fundamental (Peregrine) rogue waves 
arising at different times and spatial positions. Interesting patterns of these rogue-wave arrays are also illustrated. 



2 General rogue-wave solutions 



In this article, we consider general rogue waves in the focusing NLS equation 

iu t = u xx + 2\u\ 2 u. (1) 

Rogue waves are nonlinear waves which approach a constant background at large time and distances. Notice that Eq. 
([I]) is invariant under scalings x — >• ax, t — » a 2 t, u — > u/a for any real constant a. In addition, it is invariant under 
the Galilean transformation u(x, t) — > u(x — vt, t)exp(—ivx/2 + iv 2 t/A) for any real velocity v. Thus we only consider 
rogue waves which approach unit-amplitude background at large x and t, 

u(x, t) -> e~ 2lt , x, t -> ±oo. 

Then under the variable transform u — > ue~ 2lt , the NLS equation (TTJ becomes 

iu t = u xx + 2(|u| 2 - l)u, (2) 

where 

u{x,t) — > 1, x,t — > ±oo. (3) 



The rogue waves are described by rational solutions in the NLS equation. In order to present these solutions, let 
us introduce the so-called elementary Schur polynomials S n (x) which are defined via the generating function, 

oo / oo \ 

n=0 \k=l / 

where x = (x±, X2, ■ ■ •)• For example we have 

S (x) = 1, S 1 (x)=x 1 , S 2 (x) = \x\ + x 2 , S 3 (x) = \x\ + xix 2 + x 3 , 

2 o 

It is known that the general Schur polynomials give the complete set of homogeneous-weight algebraic solutions for 
the Kadomtsev-Petviashvili (KP) hierarchy [TBI IT?] . 
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Theorem 1. The NLS equation @ under the boundary condition © has nonsingular rational solutions 

0"! 



u = 



(4) 



where 

° n = 1<S<7V ( m 2i-l,2i-l) > ( 5 ) 

the matrix elements in a n are defined by 



mf= £ zWvff, $£> = ^Y j a k S i . v - k { X + {n) + vs), = _L W + vs), (6) 

f=0 fc=0 " i=0 

a^, (fc = 0, 1, • ■ ■) arc complex constants, and x^(n) = (xf(n), x% , • • ■), s = (si, s 2 , ■ ■ ■) are defined by 



+ / 1 ± x T 2 fc ^ 

x 1 (n) = x=F2zt±n - -, ir fc = r k , (k>2), 

°° A °° 2 A 

J] r fc A fe = ln(cosh -), ^ s fc A fc = ln(- tanh -). 



(7) 



fc=i fc=i 
The above a n can also be expressed as 

1 3 5 2N-1 



E. £.~ E .j,(«m,)j,(«m,). <»> 

where we further define 



fl=0 V2=V\ + 1 v 3 =v 2 + l v N =v N -i + l 



*l n) =0, *£>=(), (i<z/). (9) 

Before deriving these rogue wave solutions in this theorem, we give some comments. In the above definitions 
of rfe and s k , since the generators are even functions, all odd terms are zero, i.e., n = = rs = • ■ • = and 
Si = S3 = S5 = ■ • • = 0. The even-term coefficients are 

III 1 7 31 



Tl = ~, r 4 = -T7TK> r 6 — 7T7^7;> S2 = -T7T, s 4 = 7-7-77:, « 6 



192' u 2880' ' 12' 1440' 90720' 

In the solutions, a k are complex parameters. We will show in the appendix that without any loss of generality, we can 
set 

da = 1, Cl2 = CbA = • • ■ = Qeven = 0. 

In addition, by a shift of the x and i axes, we can make a\ = 0. Thus, these solutions have N — 1 irreducible complex 
parameters, a 3 , a 5 , • • •, o^jv-i- 



3 Derivation of general rogue-wave solutions 

In this section, we derive the general rogue- wave solutions given in Theorem 1. This derivation utilizes the bilinear 
method in the soliton theory [15]. The outline of this derivation is as follows. The NLS equation © is first transformed 
into the bilinear form, 

(^ + 2 )/-/ = 2M 2 , 
{Dl - iD t )g ■ f = 0, 
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by the variable transformation 

u= 9 - f (11) 

where / is a real variable and g a complex one. Here D is the Hirota's bilinear differential operator defined by 

P(D X , D y ,D t , ■ • -)F(x, y,t,---)- G(x, y,t,-- •) 

= P{d x - d x ,,d y - d yl ,d t -d t >,-- -)F(x,y,t,-- -)G(x' ,y' ',t' ', ■ ■ ■)\ x , =XtV , =y ^ =tt ..., 

where P is a polynomial of D x , D y , D t , ■ ■ ■■ Then we consider a 2+1 dimensional generalization of the above bilinear 
equation, 

(D x D y + 2)f-f = 2gh, 
(Dl - iD t )g • / = 0, 

where h is another complex variable. This is in fact the bilinear form of the Davey-Stewartson equation, which is a 
2+1 dimensional generalization of the NLS equation. We first construct a wide class of solutions for Eq. (fT2"]) in the 
form of Gram determinants. If the solutions /, g and h of Eq. (|12p further satisfy the conditions, 

(d X - dy)} = Cf, (13) 

/ : real, h = g, (14) 

where C is a constant and the overbar ~ represents complex conjugation, then these solutions also satisfy the bilinear 
NLS equation (fT0|) . Among the determinant solutions for the 2+1 dimensional system (fT2|) . we extract algebraic 
solutions satisfying the reduction condition (fl"3|) . Then such algebraic solutions satisfy both (fT2|) and (|T3|) . i.e., they 
are solutions for the 1+1 dimensional system, 

(D 2 x + 2)f-f = 2gh, 

(15 

(D 2 X - iD t )g ■ f = 0. 

Finally we impose the real and complex conjugate condition (|14[) on the algebraic solutions. Then the bilinear system 
(fl~5j) reduces to the bilinear NLS equation (fT0|) . hence Eq. (fTTj) gives the general high-order rogue wave solutions for 
the NLS equation 



The execution of the above derivation will involve some novel techniques which are uncommon in the bilinear solu- 
tion method [15] . It is known that the bilinear equations of the NLS hierarchy admit homogeneous- weight polynomial 
solutions given by the Schur polynomials associated with rectangular Young diagrams [T5]. However those solutions 
do not satisfy the complex conjugation condition h = g in general, since the Schur polynomials g and h in |18j have 
different degrees unless the Young diagram associated with / is a square. In the case of a square-shape Young diagram 
for /, h can be equal to — g (but not g) and the equation becomes the defocusing NLS equation. To construct rational 
solutions for the focusing NLS equation (|10j) . it is crucial to consider wcight-inhomogencous polynomials. In order to 
satisfy the reduction condition (flUf as well as the complex conjugate condition (|14p . we need a specific combination 
of Schur polynomials as given in Theorem 1. 



Next, we follow the above outline to derive general rogue-wave solutions to the NLS equation ([2]) under the 
boundary condition ([3]). 



3.1 Gram determinant solution for the 2+1 dimensional system 



In this subsection, we first derive the Gram determinant solution for the 2+1 dimensional bilinear equations (| 12[) . 
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Lemma 1. Let rn^\ tpr 1 ' and ifrj 1 ' be functions of x%, x-i and X-\ satisfying the following differential and difference 
relations, 

d x _ im ^ = -^r l) ^r l \ (i6) 

m (j+D =m (j) +¥ ,(n) V; («+l) ) 

^^ n) = ^" +fc) , ^Vi n) = -^ ra - fc) , (A; = 1,2,-1). 

det [m^), (17) 

l<i,j<N v 



Then the determinant, 



satisfies the bilinear equations, 

(AciAz_i — 2)t„ • T„ = — 2T n+ iT n _i, 

(18) 



Proof. We have the differential formula of determinant, 



<9 X det (dj,-) = AjJ^f 
l<i,j<N ^— ' 



• j ■ 



(19) 



and the expansion formula of bordered determinant, 

b, 



det 



Cj d 



) , A, , /),c, + d (let (ay), 



where Ay is the (i, j)-cofactor of the matrix (ay). By using these formulae repeatedly, we can verify that the 
derivatives and shifts of the r function (fl7|) are expressed by the bordered determinants as follows, 



l- in- 
to t% 

*J 

,(«) 
(n) 

(n) 
U 

(n) 



(») 



-Vj n) o 



-V 



(n+1) 


+ 


(n) 
"»« 


(n) 







, (n-1) 





(n+1) 




(n) 


(«) 













(n) (n-1) 







(n+1) 



{d Xl d x _ 1 - l)r„ 



(n) (n-1) (n) 







«7 
(n+1) 



Tn+1 



9 



-t 

(n) 





-1 



(n+1) 



5 



T n -1 

d Xl T n +i 



dljn+l 
d X2 T n+l 

From the Jacobi formula of determinants, 

bi Ci 



(n) (n-1) 



i j 

in) 



An) 
h ij 

(n+1) 



An) 
(n+1) 



(n+1) 


(n+2) 







An) 
h ij 

(n+1) 



(n+2) 
<Pi 



(n) 


(n) 


(n+1 












1 





(n) 


(n) 
(Pi 


(n+1 


-vj n) 










1 






e / 
t/j ft, fc 



X | Ojj | 



we immediately obtain the identities, 





X 




b, 






bi 


X 






9j k 


dj 


e 




9j 


h 




./' 


n &x i T n X 


d x _ 




-(- 




-i)(- 


-T„- 


u), 





^(<9^ - d^Kn+i x r n = d Xl r n+1 x d Xl r n - T n+ i^(d xi + d X2 )r n , 



which are the bilinear equations (fT8]l . This completes the proof. 



□ 



Since the matrix element m„-? is written as 



the determinant (|17[) is often called the Gram determinant solution. Let us define 

f = T , g = n, H = t-i, 

then these are the Gram determinant solution for the 2+1 dimensional system, 

{D xl D x _,-2)f- f = -2gh, 
(D 2 xl -D X2 )g-f = 0, 

which is nothing but the bilinear equations (|12p by writing x\ — x, X2 = —it and X-\ = —y. 
3.2 Algebraic solutions for the 1+1 dimensional system 

Next we derive algebraic solutions satisfying both the bilinear equations ([12")) and the reduction condition (TT3]), hence 
satisfying the 1+1 dimensional system (|15[) . These solutions are obtained by choosing the matrix elements appropri- 
ately in the Gram determinant solution in Lemma 1. 
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Lemma 2. We define matrix elements mf^ by 



p=l, g =l 



mW = 



-(— ) n e e+?7 , £=pxi+p 2 x 2 , r] = qxi-q 2 x 2 , 



p + q q 

where Aj and are differential operators with respect to p and q respectively, defined as 

Ao = ao, 

Ai = aopdp + oi, 

^2 = 'rripdp) 2 + axpdp + a 2 , 



fe=0 



(i - fc)! 



and 



-Bo — &o j 

Si = 6o?9 g + 61, 

B 2 = ^(qd q ) 2 + b iq d q + b 2 , 



and afc and 6; are constants. Then the determinant 



j-t 



r n = det 

l<i,j<N 



( (") 

I m 2i-l,2j-l 



>) 
'31 



(n) 

(«) 
7 -33 



» 



» 



"2iV-l,l "'2iV-l,3 



satisfies the bilinear equations 



(D 2 Xi + 2)t„ • t„ = 2r n+ ir„_i, 
(D^ _ D X2 )r n+1 • T n = 0. 



(n) 
TO 1,27V-1 

(n) 
TO 3,2AT-1 



An) 

l 2N-l,2N- 



Proof. First let us introduce rrS n \ tp( n ' and by 



1 -(--) n et+\ (p^=p n e^, 4> (n) = {-q)- ,l e\ 



where 



p + q q 



? 1 2 - 1 2 

£ = — ac— 1 + p x 2 , i] — -x-i + qxi - q x 2 . 



Obviously these functions satisfy the differential and difference relations 

d xl mS n) = (p {n) ^ n \ 

d Xk (p^ = (p^+ h \ dsrfW = (k = 1,2,-1). 

Therefore, by defining 

we see that these ^|™' and obey the differential and difference relations (fTS)) since the operators Aj and 

commute with differentials <9 Xfc . Lemma 1 then tells us that for an arbitrary sequence of indices 12, • • ■ , in] Ji, J2, • ■ ■ , On), 
the determinant 

f„ = dct (rhf', ) 
satisfies the bilinear equations (JTHJ)- For example, 

f " = i<S^H"-i, 2 ,-i). 

is a solution to Eq. (|18p. where p and q are arbitrary parameters. 

Next we consider the reduction condition. From the Leibniz rule, we have the operator relation, 

(pdp) m (p+-) = E ( 7 )(p+ (-i)'-)(p^) ro ~'. 



thus we get 



i-l 



/ — — I> — 



and similarly 



j 



(n) 

By using these relations, we find that fh\ ■ satisfies 



(d xi + flx-JmW - AiBj(d xl + x _JmW = AiB^p + q + - + -)m^ 

P 1 

= J2^p+ (-1)* Vi-fc^m (n) + E^+ (-i)'-)M-ift (fl) 

k=o 1 z=o y 

= E > + +E^ + (-1)' 

fc=0 ' F i=0 ' y 



<s 



Now let us take p = 1 and q = 1. Then m 



p=i,g=i 



satisfies the contiguity relation, 



fc=0 
fc:even 



£:cvcn 



p=l,q=l 



(24) 



By using the formula (|19[) and the above relation, the differential of the determinant, 



f „ = det 

l<i,j<N 



is calculated as 



iV N 

(d Xl + d x _ x )f n = A v( d *i + d *-i) 

i=i j=i 

/ 2i-l 



~ (») 



- (n) 

m 2i-l,2j-l 



N N 



A': eve ii 



JV 2i-l 



V 



,= l, g =l 
p=l,g=l 

2i-i 



p=l,g=l ' — ' l\ 
Z:cvcn 



p=l,q=l 



JV 



2 E E pE^«i-w-. 

i=i k=o ' j=i 

fc:even 



where Ay is the (i, j)-cofactor of mat I m^i-i 27 _i 



2 V - m (n) 

7] " J 2i-l,2j-l 

Z=0 
Z:evcn 

JV 2j-l TV 



j'=l Z=0 i=l 
i:even 



p=l, 9 =l 



l<i,j<N 



=1.0=1 



In the first term of the right-hand side, only the 



term with k = survives and the other terms vanish, since for k = 2,4, • • •, the summation with respect to j is a 
determinant with two identical rows. Similarly in the second term, only the term with I = remains. Thus the right 
side of the above equation becomes 



TV N 

2 E E 

i=i j=i 



- («) 



N N 



=1,0=1 



j=l i=l 



=1,9=1 



4iVf„. 



Therefore f „ satisfies the reduction condition 



(25) 



Since f„ is a special case of f„, it also satisfies the bilinear equations (|T8|) with r„ replaced by f„. From (jT8|) and ([25 
we see that f„ satisfies the 1+1 dimensional bilinear equations 

(P Xx + 2)f n ■ f n = 2f „+ir „-i, 
- £>x 2 )r„+i • f„ = 0, 



which are the same as Eq. (f2"3|). Now we can take x_i =0, then m 



- («) 



P=l,9=l 



and f n reduce to mS'f 1 and r n in Lemma 



2, and this r„ satisfies the bilinear equations (|2"5)l . This completes the proof. 



□ 



The above proof uses the technique of reduction. The reduction is a procedure to derive solutions of a lower 
dimensional system from those of a higher dimensional one. By using the reduction condition (|25[) . the derivative with 
respect to a variable x_i is replaced by the derivative with respect to another variable X\. Then in the solution, x-\ 
is just a parameter to which we can substitute any value (such as zero as we did above). 
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It is remarkable that the determinant expression of the solution (j2"2"]l has a quite unique structure: the indices 
of matrix elements, which label the degree of polynomial, have the step of 2. This comes from the requirement of 
the reduction condition, i.e., since the contiguity relation (|24p relates matrix elements with indices shifted by even 
numbers, we want such a determinant to satisfy the reduction condition. This type of Gram determinant solutions 
has not been reported in the literature to the best of the authors' knowledge. 

From Lemma 2, by writing x\ = x and X2 = —it, we find that / = To, g = T\ and h = r_i satisfy the 1+1 
dimensional system (|15[) . 



3.3 Complex conjugacy and regularity 



Now we consider the complex conjugate condition (| 14[) and the regularity (nonsingularity) of solutions. This complex 
conjugate condition now is 

ro : real, r_i = fi. 

Since x\ = x is real and x 2 = —it is pure imaginary in Lemma 2, the above condition is easily satisfied by taking the 
parameters and bk to be complex conjugate to each other, 

b k = o fc . (26) 

In fact, under the condition (f2"6) we have 



(n) (n) 

m\ 3 = m\ 3 



(-«) 
= m) t , 



and therefore 

Tn — T_ n . 



Under condition (|2"6")l , we can further show that the rational solution u = g/f = t\/tq is nonsingular, i.e., 
To is nonzero for all {x,t). To prove it, we notice that / = tq is the determinant of a Hcrmitian matrix M = 

mat ( TnS_i a—i I ■ ^ 0T an y non-zero column vector v ~ (vi, V2, ■ ■ • , vn) t and v being its complex transpose, we 

Ki,j<N V ' J J 



,3 

have 



p=l,g=l 



vMv = ^ v i m 2°i-l,2j-l V j = ViV i A -2i-iB 2 j-i—^—e i+v 

i.j = l i,j=l P 1 

N . x . x N 

= VjVj Azi-ilhj-i / e i+r *dx = I ^2 v. l v J A2i-iB 2 j-ie !i+v \ p=1 q=1 dx 

i,j=l P=l,«=l »,3=1 



p=i 



rfa; > 0, 



which proves that the Hermitian matrix M is positive definite. Therefore the denominator / = dct M > 0, so the 
solution u is nonsingular. 

It is noted that the above proofs of complex conjugate condition and regularity condition arc quite easy. This is 
an advantage of the Gram determinant expression of solutions (as compared to the Wronskian expression) . 
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Summarizing the above results, we obtain the following intermediate theorem on rogue-wave solutions in the NLS 
equation. 



Theorem 2. The NLS equation @ has the nonsingular rational solutions, 

n 



To 



where 



where the matrix elements are defined by 



l<i,j<N 



>) 

'22-1,2^-1 



l 1 

(n) a k a l 



^ (*-*)! (7-0! 
and afe are complex constants. 



p + q q 



(27) 
(28) 

(29) 



,=l, g= l 



3.4 Simplification of rogue-wave solutions 

Finally we simplify the rogue- wave solutions in Theorem 2 and derive the solution formulae given in Theorem 1. The 
generator Q of the differential operators (pd p ) k (qd q ) 1 is given as 

°° °° K k X 1 

Q = E E 'k\~n^ pdp ^ k ^ qdq ^ = ex P( K P®p + ^1 d i) = CXp(«<9inp + A<?l n g), 

A;=0 1=0 

thus for any function F(p, q), we have 

gF( P ,q)=F(e K p,e x q). (30) 

This relation can also be seen by expanding its right hand side into Taylor series of («, A) around the point (0, 0). By 
applying this relation to 



we get 



thus 



m (n) = ex P ((p + g)z - (p 2 - g 2 )**) , 

P + <7 Q 

/„2k 2 „2X2\ 



0m< n > = *(__!_.)» exp ({e K p + e x q)x - (e 2re p 2 - e 2 V)**) , 

e^p + e A q e q 



An 



,= 1, 9 =1 + ^ 



e" (K - A) exp ((e K + e A - 2)a; - (e 2K - e 2A )ii) 



= (e ._ 1)(eA _i ) exp L(k A) + (e* + e A - 2)* - (e 2 * - e 2X )it - In ^ + ^ + 1 A . 

1 (e* + l)(e A +l) V ' 

In the most right-hand side, the exponent is rewritten as 

"(« - A) + E fc[ ( x " 2fc ^) + E 7f + 2 ' ft ) - 2 ~ 2 " ln ( C ° Sh 2 C ° Sh 2 ) = ^ + E ^ 
fe=i ' i=i ' ^ ' fe=i i=i 
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where and x l are defined in ([7]), and the prefactor is rewritten as 



E 



v=0 



(e"-l)(e A -l)\" ^f K \ 



-i V(e K + l)(e A + l) 



E 

i/=0 



,4 k A 

exp [ v In ( — - tanh — tanh — 



i/=0 



= E(x) exp(^^ Sfc ( K fe + A A 



fc=i 



where is defined in ([7]). Therefore we obtain 



m (n) 



I T 

p=l,q=\ „ = q 



kA\ 

4 



exp ^( 

\fe=i 



i=l / 



and taking the coefficient of k A of both sides, we find 



1 1 
m(") fc!Z! 



(pd P T(qd q ym 



min(fe,Z) 



= ^ — Sfc-^af 1 " + Vs)Sl- u (x +VS). 



p—l t q—l V=Q 

Using the above results, the matrix element of the Gram determinant is then calculated as 

i j ux\n[i—k,j—V) 

E E a ^ 



i 



(n) J 



X] ^j5' l -fc- I y(^ + + vs)Sj-i- v (x +vs) 
E ^EE 8 ^'^'-^^ + ^s)S' i ,_i_ I ,(:E~ + i/s). 



p=i,9=i fc=o ;=o 

min (M') -, i— w J—v 



v=Q k=0 1=0 



Putting a n = T n /(m' n '| x ? _ 1 ) 7V , we obtain the determinant expression in ([5]) and ©. Finally by using §§§ and the 
formula, 



det(cijj + fcjCj) = det 

repeatedly, the determinant er„ can be rewritten into the following 37V x 3-/V determinant form, 



-Cj 1 



a n = det 

l<i,j<N 



rnin(2j-l,2j-l) 

E 



'2N-1 



*2*-1,i/*2j — l,v 



!/=0 



(") 

31 



(n) 

1,2JV-1 



(») 
10 
(n) 
11 



(») 
30 
(n) 
31 



(n) 

2JV-1,0 
(n) 

2JV-1,1 



T 2JV-1,0 V 2JV-1,1 



(n) 
3,2JV-1 



(») 

2W-1.2JV-1 



(n) 

1,2JV-1 



(n) 

3,2JV-1 



(n) 



- 2AT-1,2JV-1 

where O and / are the N x N zero matrix and 2iV x 2N unit matrix, respectively. Applying the Laplace expansion 
to the above determinant, we get 



E 

0<v 1 <u 2 <-<v N <2N-l 



(n) 
3i/i 



ll>2 



(n) 



(n) 



(n) 
3i^jv 



(n) 



(n) 



... $ 



(n) 



11/ 2 3^2 



... $ 



-2JV-l,^i ^2AT-l,!/ 2 Y 2Af-l,iyjv 

and noticing (|9|), the expanded expression (JSJ) is obtained. Theorem 1 is then proved. 



(n) 

2JV-l,i/ 2 

*1i/jv 3i/jv *2N-l,v N 
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3.5 Boundary conditions 



In order to show the boundary asymptotics ((3]), let us estimate the degree of polynomials of the denominator and 
numerator in The elementary Schur polynomial Sk(x) has the form Sk(x) = (xi) k /kl + (lower degree terms), 
where x = (xi, X2, • ■ ■)■ Thus the degree of the polynomial Sk(x +vs) in (x, t) is k and its leading term appears in the 



monomial (xf) k /k\, i.e., the leading term is given by (x^f2it) k /k\. Therefore the degrees of and are both j—v, 
and their leading terms are ao(x — 2it)j~ v j (j — v)\2 v and do(x + 2it) J ~ ly / (j — v)\2 u , respectively. Therefore both of the 
degrees of determinants dct ( 3> 2 i-i v I an d det ( *2i-i v J are gi ven by 1 + 3H h(2iV— 1) — i^ — ^ 2 fjv, 

l<i,j<N \ ' 3 J l<i,j<N V ' 3 J 

and in the expression ([ 
summation. For ja- = j 



|, the highest degree term comes from the term of V\ = 0, v-i 
- 1, we have 



I'N 



N - 1 in the 



det 

l<i,j<iV 



(n) 

2i-l,3-l 



ap(4) 3 
3! 



ao(x+) 2JV - 
(2N-1)! 



g%(x+) N ( N+i y 2 
1!3! • • • (2N - i)!2^(JV-i)/2 

1 1 
13 3-2 



1 2JV-1 (2JV-l)(2JV-2) (2JV 
+ (lower degree terms). 

The above determinant is equal to 



det 

l<i,j<N 



"11 
2 

1 + \2 

a {xj) 







On 
2^ 


2!2 


2^ 

ao(x+) 2N - 3 


(2iV-2)!2 


(2JV-3)!2 2 


(27V-4)!2 3 












3-2-1 






ao(x£) 
7VI2™- 1 



(lower degree terms) 



l)(2AT-2)(2JV-3) 



(2AT-l)(2JV-2)---(AT + l) 




= det ((2* - IV' -1 ) , 



which is the Vandermonde determinant. Thus we obtain 

= S'.'.S-i ) )! ffl °^ ) ^ 1>/2 + (iowcr dcgrce tcrms) ' 



det 

l<i,j<N 



and similarly 



det 

l<i,j<N 



(2N- 1)! 

^w^ ~ a ^^ )N(N+1)/2 + (lower degree terms) - 



Consequently, the leading term of a n is given by 

/om.-^-in 2 

Vl!3!---(27V-l)!y 1 v ; 
which is independent of n. Hence u = ci/ao satisfies the boundary condition J3]). 



4 Solution dynamics 

In this section, we discuss the dynamics of these general rogue wave solutions. 
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To obtain the first-order rogue wave, we set N = 1 in Theorem 1. In this case, 



(0, 



(x - lit - - + a x ){x + 2it - i + oi) + 



1 3 1 

m\{ = (x - 2it + - + a{){x + 2U-- + a x ) + -, 



hence the first-order rogue wave is 



u(x, t) 



(a; - 2rf + | + oi ) (a; + 2it - § + a x ) + i 
(a: - 2ii - ■§ + ai)(x + 2ii - | + a x ) + ± 



(31) 



Clearly, the complex parameter ai in this solution can be normalized to zero by a shift of x and t, as we have mentioned 
before. After setting a\ = 0, this first-order rogue wave can be rewritten as 



u(x, t) = 1 — 



4(1 - Ait) 



1 + Ax' 



16t 2 ' 



(32) 



where x = x — 1/2. This rogue wave was first obtained by Peregrine [7], see also jS). Its maximum peak amplitude is 
equal to 3, i.e., three times the background amplitude. 

To obtain the second-order rogue waves, we take N = 2. In this case, 

(i) (i) 
m i3 



;i) (i) 

m 31 m 33 



(o) W) 

(0) (0) 

™>3i rn 3 3 



(33) 



From the previous discussions, we will set a\ = a-i = 0. Then the general second-order rogue wave can be obtained 
from (1531 as 



u = 1 + 



1>' 



(34) 



where 



ix 2 t 2 - 160t 4 ) 



24{(3x - 6x 2 + Ax 3 - 2x 4 - 48t 2 + A8xt 2 
+«*(— 12 + 12.x - 16x 3 + 8a- 4 + 32t 2 - 64a;t 2 + 64x 2 t 2 + 128t 4 ) 
+6a 3 (l - 2x + x 2 - 4ii + 4ia;t - At 2 ) + 6a 3 (-x 2 + Aixt + At 2 )}, 



4) = (9 - 36a; + 72a; 2 - 72a- 3 + 72a; 4 - 48a; 5 + 16a; 6 ) 

+96t 2 (3 + 3a; - 4a; 3 + 2a; 4 ) + 384^(5 - 2a; + 2a; 2 ) + 1024i 6 

+24(a 3 + a 3 )(3x 2 - 2x 3 - 12t 2 + 2Axi 2 ) + 48i(a 3 - a 3 )(3i + 6xt - 6x 2 t + St 3 ) + 144a 3 a 3 , 

and a 3 is a free complex parameter. We have found that the maximum of |u(a;,i,a 3 ) is equal to 5, and it is obtained 
when 

-1/12. 



At this a 3 value, the solution is 



u m (x,t) = 1 + 



(35) 
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where 



<t> m = 9 - 72i 2 - 48i 4 - 864t 2 - 3840t 4 - llb2x 2 t 2 + it(-180 - 288i 2 + 192x 4 + 384i 2 + 3072t 4 + 1536zV), 

V> m = j + 27i 2 + 12i 4 + 16x 6 + 396t 2 + 1728i 4 + 1024< 6 - 288i 2 t 2 + 768£ 2 t 4 + 192x 4 t 2 , 

and x — x — 0.5. This solution is displayed in Fig. 1(a). It is easy to see that this solution is the special second-order 
rogue wave obtained by Akhmcdiev et al. [8] (after a shift in x) . Thus the special second-order rogue wave obtained by 
Akhmediev et al. is the one with the highest peak amplitude among all second-order rogue waves. At other 03 values, 
however, we can obtain rogue waves which have very different solution dynamics from that in Fig. 1(a). For instance, 
rogue waves at 03 = 5/3, —hi/2 and hi/2 are displayed in Fig. l(b,c,d) respectively. In each of these solutions, three 
intensity humps appear at different times and/or space, and each intensity hump is roughly a first-order (Peregrine) 
rogue wave (f32|) . Specifically, in Fig. 1(b), the solution features double temporal bumps (elevations) at x w —0.5 and 
a single temporal bump at x w 2.2. In Fig. 1(c), the solution first rises up and reaches a peak at (x,t) w (0.5, —0.7). 
Afterwards, the solution temporally decays at x « 0.5, but two new bumps rise at the two sides. In Fig. 1(d), the 
solution is similar to that in Fig. 1(c) but with a time reversal. Obviously, the rogue- wave dynamics in Fig. l(b-d) are 
quite different from the one in Fig. 1(a). The solution dynamics in Fig. l(b-d) resemble those reported in [9"1 [Till [Tl"] . 



(a) (b) 

/ smmk 




Figure 1: Second-order rogue waves with parameters: (a) 03 = —1/12; (b) 03 = 5/3; (c) 03 = — hi/2; (d) 03 = hi/2. 

Next we examine third-order rogue waves. In this case, we set a\ = a-i — 04 = without loss of generality If one 
takes 

03 = -1/12, 05 = -1/240, 

then the corresponding solution u m (x, t) is equal to the third-order rogue wave obtained by Akhmediev et al. [8] 
except a shift in x. This solution is displayed in Fig. 2(a). The maximum amplitude of this solution is equal to 7, 
which occurs at (a;, t) = (1/2, 0). We have found that this special rogue-wave solution u m (x, t) is also the one with the 
highest peak amplitude among all third-order rogue waves u(x, t; 03, 05). But if we take other (03,05) values, rogue 
waves with dynamics different from Fig. 2(a) will be obtained. Three of such solutions, with (03,05) = (25/3,0), 
(— 25i/3, 0) and (0, 50i/3), arc displayed in Fig. 2(b,c,d) respectively. These solutions feature six intensity humps 
which appear at different times and/or space, and each intensity hump is roughly a first-order rogue wave (|32[) . In 
Fig. 2(b), the solution exhibits triple temporal bumps at 1 « —2. double temporal bumps at x « 2, and a single 
temporal bump at x 6. In Fig. 2(c), the solution develops a single hump first. Then this hump decays, but two 
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new humps rise simultaneously at the two sides. Then these two humps decay, but three additional humps develop 
simultaneously. In Fig. 2(d), two intensity humps rise simultaneously at different spatial locations first. After they 
decay, additional four intensity humps arise at different locations and times. A remarkable feature in the rogue waves 
in Fig. 2(b-d) is the high regularity of their spatiotcmporal patterns. For instance, the pattern in Fig. 2(c) is a highly 
symmetric triangle, while the one in Fig. 2(d) is like a pentagon. These spatiotemporal patterns of rogue waves are 
different from the ones reported in [TUJ [TT] . 




-2-10123 
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• •< 
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. . . . 



-4 4 




-4 4 



Figure 2: Third-order rogue waves with parameters (03,05) as: (a) (—1/12,-1/240); (b) (25/3,0); (c) (— 25z/3, 0); 
(d) (0,5(K/3). 

The results shown above apparently can be extended to fourth- and higher-order rogue waves. By special choices 
of the free parameters (03,05,07, • • •), we can reproduce the rogue waves obtained in [8] as special cases. But other 
choices of those parameters can yield even richer spatiotemporal patterns, such as triangular patterns like Fig. 2(c) 
but with more intensity humps such as 10, 15, and so on. 



5 Summary and discussion 



In this paper, we derived general iV-th order rogue waves in the NLS equation by the bilinear method. These solutions 
were obtained from Gram determinant solutions of bilinear equations through dimension reduction and then further 
simplified to a very explicit form. We showed that these general rogue waves contain N — 1 free irreducible complex 
parameters. By different choices of these free parameters, we obtained rogue waves with novel spatiotcmporal patterns. 
These new spatiotcmporal patterns reveal the rich dynamics in rogue- wave solutions and deepen our understanding of 
the rogue- wave phenomena. We also showed that the rogue waves reported in [8] are special solutions with the highest 
peak amplitude among all rogue waves of the same order. 

We would like to point out that the new spatiotemporal patterns of rogue waves obtained in this paper may also find 
applications in other branches of applied mathematics and physics. For instance, the triangular rogue-wave patterns 
in Figs. 1(c), 2(c) and their higher-order extensions (with more intensity humps) closely resemble the spike pattern 
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which forms after the point of gradient catastrophe in the semiclassical (zero-dispersion) limit of the NLS equation 
(see Fig. 1 in [19]). The connection between these exact rogue- wave solutions and the semiclassical-NLS patterns is 
an interesting question which lies outside the scope of the present article. 
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Appendix 

In this appendix, we determine the number of free parameters in the rogue-wave solutions obtained in this paper. 
Since the solutions in Theorem 1 are derived from the ones in Lemma 2, we will examine solutions in Lemma 2 below. 

First we can factor out eto from Ai and bo from Bj . These factors cancel out in the formula (|lll) , thus we will set 
ao = bo = 1 without loss of any generality. 

Secondly, let us consider the effect of a constant shifting of (£1,2:2). By the shifting (xx,X2) — > (xi + a,X2 + /?), 
m' n ' in (f2"Tj) gets an exponential factor, 

m (n) m («) e e , 0=(p + q)a + (p 2 -q 2 )(3, (36) 

consequently the (i, j)-component AiBjm^ 1 ' in (I2TH) is also modified. Below we show that AiBjin^ is modified as 

A.Bjm^ A t B 3 (m {n) e e ) = e e A t B 3 m {n \ (37) 



where 



and 



^Et^W^ ^-Et^W'"', (38) 

V" c 1 + \ + ( 2flP Q + V^ pa + 2 k p 2 (3 \ 
u V^; c I -\ - ( o 2„ 9«- 4 <? 2 /3 qa-2 k q 2 f3 \ 
To prove (|3"T)) . we notice that for the generator Q of the differential operators (pd p ) k defined by 

G = E T\ (P 9 ^ = eX P( A ^f) = CX P( A9 lnp), (41) 
fe=0 

the relation 

gF(p,q)=F(e x p,q) 
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holds for any function F(p,q). This relation is a special case of the previous relation ([5U]) . Thus, 

e~ e G(e e F) = exp ((e A - l)pa + (e 2A - l)p 2 p) QF = exp ^ — (pa + 2 k p 2 (3) QF, 
whose coefficient of order X k gives 



Similarly we have 



i/=0 v ' 

h q d q )\e e F) = e e j2Su( X - ) 1T ^(qd q ) l ^F 



Therefore, 



A ? i? J (m(")e e )=^^a fe 6 i -^(pa p r- fc -^(^r i (m(")e e ) 

= * e EE E ^(^ ...^^. W-^ E s,(x ) J_ ( qdq y-^m^ 

k=o 1=0 m=o ^ v=o u 



which proves Eq. fl37p . 

Now we take p = g = 1. Then from Eqs. ([39 |) -(|30 ]l , we get 

&0 = o-o = 1, 01 = 01 + a + 2/3, 6o = 6o = l, 61 = bi + a - 2/3, 

Thus by a shifting of (xi, X2) — > (x\ + a, x 2 + /3) with a = — (ai + 6i)/2 and /3 = — (ai — &i)/4, we obtain di = 61 = 0. 
When this shifting is combined with shifts of higher coefficients a 2 — > a 2 , 03 — > 03, • • •, 62 => 63 - ^ 63, • • •, the 
solution r n depends on parameters (a 2l S3, • ■ • ; 62, 63, ■ ■ ■) only. In other words, by a shift of (x\, x 2 ), we can normalize 
ai = b\ = 0. 



Thirdly, from the expressions of to™- in (|20j) and the expressions of A,; and we see that in the determinant 
formula for r„ in (|22[) , when we subtract the product of the first row and 02 from the second row, and subtract the 
product of the second row and a 2 from the third row, . . . , and subtract the product of the ith row and a 2 from the 
(i + l)-th row, and then subtract the product of the first column and 62 from the second column, and subtract the 
product of the second column and 62 from the third column, etc., we can remove the parameter a 2 and 62 from the 
solution formula (|22j) . By similar treatments, we can remove all other even coefficients a 4 , ae, ■ ■ ■ and 64, be, ■ • • as well. 
In other words, we can set a 2 = 04 = ag = • • • = and b 2 = = bg = ■ ■ ■ = without any loss of generality. 

By summarizing the above results, we see that without any loss of generality, we can set 

do = bo = 1, a 2 — 04 = a 6 = • ■ ■ = b 2 = b 4 = 6 6 = • • • = 0. 

In addition, by a shift of (xi,x 2 ), we can normalize a\ = 61 = 0. Combined with the complex conjugacy condition 
bk = Sfc in (|26[) . we then find that the TV-th order rogue- wave solutions in Theorem 1 have N — 1 free irreducible 
complex parameters, 03, 05, ... , (i2N-i- 
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